Fluctuations of a dipolar condensate: a signature for rotons 
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We demonstrate that measurements of atom number fluctuations in a trapped dipolar condensate can reveal 
the presence of the elusive roton excitation. The key signature is a super-Possionian peak in the fluctuations 
as the size of the measurement cell is varied, with the maximum occurring when the size is comparable to 
the roton wavelength. The magnitude of this roton feature is enhanced with temperature. The variation in 
fluctuations across the condensate demonstrates that the roton excitations are effectively confined to propagate 
in the densest central region, realizing a density trapped roton gas. While our main results are based on full 
numerical solutions of the meanfield equations, we also develop and validate a simple local density theory. 
Finally, we consider fluctuations measured within a washer-shaped cell which filters out the contribution of 
modes with non-zero angular momentum and provides a signal sensitive to individual roton modes. 
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PACS numbers: 67.85-d,67.85.Bc 

Introduction Bose-Einstein condensates (BECs) have been 
created with the highly magnetic atoms Chromium |[T|, Dys- 
prosium |2| and Erbium [3]. The key new feature of these 
systems is a long ranged and anisotropic dipole-dipole inter- 
action (DDI), which gives rise to a rich array of new physics 
for ultra-cold gases |4|. An important prediction made in 2003 
was that a rotonic-like excitation would emerge in a conden- 
sate that is tightly confined along the direction that the dipoles 
are polarized |5| (also see [6J). To date there has been no 
experimental evidence for the existence of these excitations, 
and increasing attention is turning to developing signatures 
for detecting their presence, such as structure formation in the 
ground state density profiles |6-9J, anisotropic superfluidity 
[[TOllITl and sensitivity to perturbations fTZ'-^. 

The rotons within dipolar BECs arise when tight axial con- 
finement, of length scale Iz, causes the effective interaction in 
the weakly confined (in-plane) directions to become momen- 
tum dependent |15|. For in-plane momenta hkp < h/lz the 
DDI is repulsive, while for hkp > h/lz it crosses over to being 
attractive, energetically softening these modes and causing the 
roton to emerge with a wavelength Ajot ^ Iz- These rotons are 
sensitive to the condensate density and, due to the radial trap- 
ping, they are effectively confined as a roton gas in the high 
density central part of the condensate lfT6l . 

Here we investigate how rotons manifest in the fluctuations 
of a dipolar BEC. This is motivated by significant develop- 
ments made in BEC imaging systems to allow precise in situ 
measurements of density (e.g. see Iil71 ). Such measurements 
actually determine the atom number N^r within a finite sized 
cell a, with the cell parameters determined by the imaging 
system ifTSll . In equilibrium the number fluctuations 



5Nl = {ND - Nl 



(1) 



about the mean number N^ = {No) are crucially dependent 
on the cell size fWj : for cells smaller than all relevant length 
scales (i.e. healing length and thermal de Broglie wavelength) 
the fluctuations are Poissonian 5N'^ ~ N„, whereas for large 
cells the thermodynamic result holds SN^ ~ Na-ksTriKT, 
with n the density and kt the isothermal compressibility. 
DDIs introduce two important features for fluctuations: (1) 
the roton wavelength introduces a new length scale which is 



revealed as the cell size changes; (2) the attractive nature of 
the interaction for short wavelengths induces enhanced fluc- 
tuations, this is in contrast to the suppression of fluctuations 
expected in a condensate with short ranged s-wave interac- 
tions Ho). 
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FIG. 1. (color online) Schematic geometry of the system we con- 
sider: a radially symmetric pancake condensate realized by tight con- 
finement along z; cylinder and washer-shaped cells in which number 
fluctuations are measured. The cylindrical cell is parameterized by 
the radius R from trap center to the cylinder axis, and its diameter 
D. The washer-shaped cell is centered on the trap center, and is pa- 
rameterized by the mid-radius R and width W. For both cells the z 
extent is taken to be larger than the condensate thickness. 



Formalism The geometry of the tightly confined (pancake) 
dipolar condensate and the types of cells in which we measure 
number fluctuations are shown in Fig.fT] The two cell shapes 
we consider are cylindrical and washer-shaped cells with the 
z extent of these cells being larger than the condensate thick- 
ness in this direction. The cylindrical cell can be realized ex- 
perimentally using absorption imaging (column density along 
z) with the minimum diameter (D) being set by the spot size 
of the imaging system (e.g. see |18] |2T] |22l)- These mea- 
surements are made using a charged coupled device (CCD) in 
which the pixel size is typically much smaller than the spot 
size. Cells with more general shapes in the xy-plane can be 
realized by collecting signal from a set of pixels, e.g. larger di- 
ameter cylinders and the washer-shaped case that we consider 
later. 



The cell number operator is Ncr = f (Pr ^^ (r)^(r), with 
^(r) the bosonic field operator and a denoting the cell vol- 
ume. Within a meanfield approach, valid for weakly interact- 
ing condensates, the field operator takes the form 



*(r)« Vo(r) + Xl [^'^W"* 



■{r)al 



(2) 



In the above expression ipQ is the condensate wavefunction 
(normalized to the condensate number A^o). which satisfies 
the non-local dipolar Gross-Pitaevskii equation (GPE) |23| 



mV'o 



2m 



Vir)+ / dr'f/(r-r')|Vo(r')l' 



i^o, (3) 



where /i is the chemical potential. For dipoles polarized along 
z the inter-atomic interaction potential is of the form 



C/(r) = gSir) + 



Sgdd 1 — 3 cos^ 9 
^TT |r|3 ' 



(4) 



where the short range interaction is characterized by the con- 
tact parameter g = Anagh^ /ra, with a^ the s-wave scattering 
length. The DDI parameter is g^d = A^oZ-'m/^i with /i,„ is the 
magnetic dipole moment, and 9 is the angle between r and 
the z axis. The atoms are taken to be harmonically confined 
by a potential V{r) = ^■muj'i{p^ + A^z^) with A = w^/wp. 
The amplitudes {ui, Vi}, with energy eigenvalues ej, describe 
the linearized quasiparticle excitations, and are obtained by 
solving non-local Bogoliubov de Gennes equations 11241 . The 
quasiparticle operators satisfy [ai, a A = 6ij. 

The mean cell number is N^- = J (Prn{r), where n(r) — 
|V'o(r)P + n{r) is the total density, with 

n(r) = J2 K>^WI' + K- + 1)I«.-WI'] , (5) 



the depletion (non-condensate) density, fij = (e^'^J — l) 
and /3 = l/ksT the inverse temperature. Neglecting terms 
greater than second order in the quasiparticles, we obtain 

5Nl =J2 I ^^1 / rfr2<5n,(ri)K(r2)coth(i/3e,) , (6) 



250 

200 

,150 

'ioo 

50 h 



14 
12 
10 

3" 8 

^ 6 

4 



(a) 






^^^ 




unstable 


^^ 


^^ 


/ 


^' 






y" a 


p 


7 


\ 






stable 




•^ 







-10 



10 

Os/Oo 



20 



30 



(b) 




'■ — — r ' ' 1 ' '. ^ 


_2 




—- 


1 


A, \ 1 


2 4 6 8 


— 






^^Sss 





-^ 




1 


- ^ 


J- ^^^ ^ 


i 




T^^' 


i 


,/ 


J 






/ 


a 


13 


7 



10 



s/a-o 



15 



20 



25 



FIG. 2. (color online) (a) Stability phase diagram in interaction pa- 
rameter space for A^'o — 2.5 x 10'' '"'^Dy atoms in a A = 20 trap 
with ojp = 27r X 10.8 s~^ . The shaded region indicates where a stable 
solution exists and the gray portions of the boundary mark the bound- 
ary of density oscillation islands 16}. Crosses indicate the three sets 
of parameters we examine in this paper with Os/ao = 0, 9.04, 18.1 
as the cases we label as q, /3, 7, with ao the Bohr radius. The dipole 
length is add = rag^d/^T^f? ~ 132.6 ao for ''^''Dy. (b) Quasiparti- 
cle energies of m^ — modes as a^ is varied along the dashed path 
in (a). Vertical lines mark the cases a, /3 and 7, and the crossed cir- 
cle indicates the lowest energy roton mode.The contours of the den- 
sity fluctuation Srii for the roton shown in inset. Contour spacing is 
3 X lO^'cm"^, solid lines (dot-dashed) represent positive (negative) 
contours, while nodes are shown as thick black lines. 



with 6ni = ■00 {ui — Vi) the density fluctuation amplitude of 
the i-th quasiparticle, where we have taken i/jq to be real. 

In Fig. l2la) we show the numerically determined stability 
diagram for Nq = 2.5 x 10^ ^^''Dy atoms in a trap of as- 
pect ratio A = 20, indicating the range of contact and dipolar 
interactions where all the quasiparticle energies are real and 
positive. In the upper part of the stable region, where the DDI 
dominates, roton excitations tend to develop (e.g. see SlAJ ). 
Indeed, in Fig. |2|b) we show the energy spectrum along the 
dashed path in Fig.l2|a) where we vary Ug, as can be done us- 
ing a Feshbach resonance. As a^ becomes smaller we observe 
a radially excited roton mode soften, leading to instability for 
fls « — 3ao. We note that there is a range of roton modes soft- 
ening, including those with non-zero values of angular mo- 
mentum projection m^ (which are not shown in the spectrum) 



although the modes that soften all have a similar wavelength, 
as indicated in the inset to Fig.[2|b). In Fig.|2]we indicate three 
parameters sets that we will use here and have labeled a, /3, 7 
(differing by the value of a^). The purely dipolar case a has 
the most prominent rotons, while cases (3 and 7 (progressively 
further from the stability boundary) may be more convenient 
to explore experimentally. 

Cylindrical Cells We present results for 6N^ for the purely 
dipolar case a in Fig. l3] These results are calculated using 
Eq. (J6]l in normally ordered form and include contributions 
from ^ 10'* quasiparticle modes. We emphasize that our cal- 
culations are fully three-dimensional and we have performed 
tests with up to 10^ modes finding changes in SN^ of <0.7%. 
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FIG. 3. (color online) SN^ as a function of cylinder position R for 
various diameters D for case a. (a.)T — with a depletion of iV = 
171. (b) T = 5.2 nK with iV = 1570. Insets: Condensate, depletion 
and fluctuation density. Densities evaluated at z = and np is 
calculated using modes with £_, < 2Ghujp. 



The total depletion is TV = / d^r n « 1570 (i.e. 6% of Nq) 
at the highest T considered. The peak in the depletion den- 
sity [see inset to Fig. [3jb)] near trap center reveals the pres- 
ence of the confined rotons modes. To emphasize that the 
rotons have a larger effect on the fluctuations than is sug- 
gested by the depletion, we define and calculate ?l fluctuation 
density by the coiTelator np (r) = ■\/G'(^)(r) — n(r)2, where 
G(2) ^ ^(•^t>)2(-^^2N) (c.f. quasicondensate definition in ESl). 
While np is generally less than n in a repulsive condensate 
with contact interactions, for the dipolar gas it is significantly 
enhanced [see inset to Fig.lSlb)] |26|. 

At T = [Fig. [Sja)] and for i? = (cells at trap center) 
the number fluctuations are almost Poissonian (6N^/Nc « 1) 
for the smallest cell diameter D — 0.239 /im, increase to be- 
come super-Poissonian for 13 ~ 2 /xm, and then decrease to 
sub-Poissonian for D > 3 fim. The dominant contribution to 
the fluctuations in Eq. (J6]) arises from modes with wavelength 
comparable to the cell width 1191 . The dominant modes for 
the smallest size cells are the high energy (free particle-like) 
excitations, while for the largest cells it is the long wavelength 
repulsively interacting phonon-like modes that suppress fluc- 
tuations. For intermediate cell sizes comparable to half the 
roton wavelength Arot/2 ~ 1.5 /im [see inset to Fig.l2tb)] the 
fluctuations are enhanced by the effective attractive character 
of the DDI for modes at this scale. This roton peak in 5N^ 



is a key feature of dipolar condensates and can be used as a 
signature of the emergence of roton excitations. 

Temperature plays an important role in BEC fluchiations, 
and the results in Fig. IStb) are for a case typical of the low- 
est temperatures obtained in experiments (« 13% of the criti- 
cal temperature) [27|. These results show that thermal effects 
wash out the sub-Poissonian fluctuations in large cells, while 
the roton peak in 6N^ (occurring for D ^ Ai.ot/2) is signif- 
icantly enhanced. For both temperature regimes in Fig. |3]the 
height of the roton peak decreases with R as the cell is moved 
away from the central region, because the rotons are effec- 
tively confinement to trap center 1161 . 
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FIG. 4. (color online) Peak fluctuations (at R — 0) within a cylindri- 
cal cell (a) as a function of temperature at diameter D = 1.67 /^m, 
and (b) as a function of cylinder diameter at T = 5.2 nK. Each 
graph has results for three different values of contact interaction (see 
Fig.[2](. Solid: full calculations and dashed: LDA result (see text). 



Figure Wa) examines the fluctuations at i? = for a cylin- 
der of fixed diameter D = 1.67 /im, chosen to be sensitive to 
the roton excitations, as a function of temperature and for sev- 
eral values of the contact interaction. In all cases the plateau at 
low T represents the quantum fluctuations associated with the 
rotons, while the transition to a rising line coincides with the 
thermal activation of these modes, occurring when fc^T is of 
the order of the lowest roton energy. With reference to spectra 
in Fig. l2tb) we note that the energy of the lowest roton in- 
creases as Qs increases. In Fig.|4|b) we show the dependence 
of SN^ on cell diameter. Importantly, the results with larger 
Os (where the roton energy is higher), and the roton peak in 
6N^ at D ^ ^Arot is smaller, but still discernible. 

In Fig. Wb) we also present results of a simplified theory 
based on the local density approximation (LDA) for calculat- 
ing the fluctuations in the quasi-two-dimensional (quasi-2D) 
regime, valid when /c^T < huj^- This LDA theory makes use 
of the GPE wavefunction but avoids the computational chal- 
lenging step of calculating all the quasiparticles and integrat- 
ing them over the cell region a, as required to use Eq. (|6]l. We 
first note that for a uniform quasi-2D system the fluctuations 
in a disk of radius D is 



6N' 



n2D d k 



JK\Dkp) 



^(kp), 



(7) 



with S'(kp) the static structure factor, n2D the areal density, 
and the Bessel function term arising as a geometric factor for 
the cell [19J. 
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FIG. 5. (color online) Fluctuations versus washer radius R and for 
various washer widths W for the a^ =0. (a) T = results. Inset: 
5ni for the lowest roton mode [repeated from inset of Fig.|2jb)]. 
(b) T — 5.2 nK results. Inset: The dashed curve shows the number 
fluctuations calculated using only the lowest roton mode, compared 
to the full calculation (solid). 



The basis of our approach for extending this result to 
the trapped system is to use i/jq to calculate a local Bo- 
goliubov dispersion relation eB(kp,p) at each radial posi- 
tion p = {x,y) (see [14| for details), from which a lo- 
cal static structure factor can be obtained as ^(kpjp) = 



coth (I/^eb) . The static structure factor for the sys- 
tem within the cell cr can then be computed by the density 
weighted LDA sum Sa(kp) — J d?pn2D{p)S(kp, p), where 
n2D{p) = J dz \ipo{p, z)\'^. The number fluctuations for the 
trapped system is then given by Eq. (|7]l with the replacement 
5(kp) — > S'cr(kp). The comparison in Fig.|4tb) demonstrates 
that the LDA results provide a useful quantitative estimate. 
Washer-shaped Cells It is interesting to consider more gen- 
eral cells shapes for measuring SN^. Here we examine the 
non-simply connected washer-shaped cell (see Fig. [Til. We 
show 6N^ for this cell versus washer radius R and for vari- 
ous washer widths at both zero [Fig. |5ja)] and finite temper- 
ature [Fig. |5jb)]. The most striking feature of these results 
is the emergence of fine oscillations in SN^ as R is varied. 
This arises because the washer geometry only allows rriz = 
modes to contribute [28 1, thus only a small number of modes 
contribute to the result. Indeed, in the inset to Fig. l5lb) we 
compare the full result for SN^ to that calculated by restrict- 
ing the sum in Eq. (J6]) to just the lowest energy roton mode, 
and find reasonable agreement, with the difference being due 
to other excited m^ = modes. The broad envelope of these 
oscillations is similar (as a function of R) to that observed for 
the cylindrical cells in Fig. [3] (reflecting roton confinement). 
Conclusions In this paper we have examined the number 
fluctuations that occur within small cells of a dipolar BEG. 
We have shown that these fluctuations provide a large signal 
indicating the presence of roton excitations, even when the de- 
pleted density associated with these excitations is small. Our 
results show that the fluctuations are highly dependent on the 
cell size relative to the roton wavelength, and should be prac- 
tical to measure in experiments with imaging that can resolve 
cells of size ^ 2 /im. This resolution has been approached 
in similar style quasi-2D experiments with Cs [18] , and well- 
exceeded with other quantum gas microscope (e.g. ||29| ) and 
in scanning electron microscopy |'30l experiments. While we 
have focused on a physically realizable regime with Dy, we 
emphasize that Cr and Er systems are also suited to this type 
of measurement. 
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